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ABSTRACT 

Epigenetics are the heritable changes in gene ex- 
pression or cellular phenotype caused by mecha- 
nisms other than changes in the underlying DNA se- 
quence. After mitosis, it is thought that bookmark- 
ing transcription factors remain at promoters, reg- 
ulating which genes become active and which re- 
main silent. Herein, we demonstrate that poly(ADP- 
ribose)polymerase-1 (PARP-1) is a genome-wide epi- 
genetic memory mark in mitotic chromatin, and we 
further show that the presence of PARP-1 is abso- 
lutely crucial for reactivation of transcription after 
mitosis. Based on these findings, a novel molecular 
model of epigenetic memory transmission through 
the cell cycle is proposed. 

INTRODUCTION 

Mitotic chromatin is transcriptionally inactive (1). The ma- 
jority of transcription factors and ribonucleicacid (RNA) 
polymerases are excluded from mitotic chromatin, albeit 
a few are retained (2-7). Following the completion of cell 
division, regulatory proteins need to re-engage their ge- 
nomic targets to restore appropriate gene transcription 
states (8-12). Epigenetic marking of mitotic chromatin is 
involved in precise post-mitotic restoration of proper tran- 
scriptional patterns to avoid disastrous regulatory conse- 
quences (13,14). Such marks include certain histone and 
DNA modifications that remain in mitotic chromatin (15- 
18). Deoxyribonucleic acid (DNA) methylation is known to 
suppress transcription upon completion of mitosis (19,20), 
whereas relation between specific histone modifications and 
gene expression states remains somewhat illusive. To reset 
different active states of genes after mitosis, transcription 
factors and their attendant co-regulators must find their ap- 
propriate sites in transcriptionally silent chromatin. Mech- 
anisms responsible for re-engagement of transcription fac- 
tors and their co-regulators in chromatin after cell divi- 
sion are in many respects similar to mechanisms respon- 
sible for chromatin programming during cell differentia- 
tion, when regulatory factors engage transcriptionally silent 



genes in progenitor cells and activate expression de novo. 
Histone modifications can help recruit transcription factors 
and co-regulators to target sites in chromatin (21-23). Early 
embryonic cells already possess distinct histone modifica- 
tion patterns at a subset of silent genes, apparently mark- 
ing genes cinched to be activated (24). Nevertheless, these 
known modifications of histones and DNAs likely represent 
only a subset of epigenetic marks (25). 

In a living cell, the activation of chromatin-associated 
PARP-1 leads to assembly of poly(ADP-ribose) (pADPr) 
via transfer of ADP-ribose to various substrates, including 
existing pADPr chain, PARP-1 itself and other acceptor 
proteins. Subsequent accumulation of pADPr neutralizes 
positively charged groups, mediating chromatin deconden- 
sation and stimulating transcription (26-30). pADPr chain 
averages 80 moieties, although length variation is substan- 
tial (31). PARP-1 is active as a dimer (32,33). A cluster of 
10-28 Glu residues located near the center of the PARP-1 
dimerization domain serves as the major acceptor of ADPr 
moieties. Auto-ADPr-modification of PARP-1 affects its 
dimerization domain, causing PARP- 1 dimers to dissociate 
from each other and from active chromatin and lose activity. 
This negative feedback loop mediated by automodification 
limits the time during which PARP- 1 remains active. 

In steady state conditions, PARP-1 is unevenly dis- 
tributed in the chromatin and often found near the pro- 
moter regions of specific genes (34,35). We hypothesized 
that PARP-1 is an epigenetic memory factor, capable of 
marking specific chromatin regions during mitosis to be re- 
activated at the start of interphase. We report the results 
of whole-genomic comparative analysis of PARP- 1 distri- 
bution in untreated asynchronous human embryonic kid- 
ney (HEK293) cells and cells synchronized in the mitotic 
phase by nocodazole treatment. This experimental system 
allows accumulation of up to 99% of a total cell population 
in mitosis (36). These cells remain alive, and upon noco- 
dazole removal immediately mitosis (Supplementary Fig- 
ure S1A-D). Our experiments demonstrated that PARP-1 
establishes stable epigenetic marks in metaphase chromatin 
at the transcription start sites of a large subset of genes and 
that the presence of these marks is necessary for transcrip- 
tional restart after mitosis. 
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MATERIALS AND METHODS 

Cell lines and cell culture 

Human Embryonic Kidney (HEK293) cells (catalog no. 
CRL-1573) were purchased from American Type Culture 
Collection (Boston, MA, USA) and cultured in Dulbecco's 
modified Eagle's medium (DMEM) supplemented with an- 
tibiotics (streptomycin, 100 mg/ml; penicillin, 100 IU/ml), 
and 10% (w/v) newborn calf serum (NCS) at 37°C in a 
humidified atmosphere containing 5% CO2. Nocodazole, 
olabarib and Epidermal Growth Factor (EGF) were used to 
treat cells in specific experiments. Short hairpin (sh) RNA 
control and shRNA PARP-l-transfected cells were syn- 
chronized using nocodazole and released in the presence 
or absence of olaparib. PARP- 1 protein in knockdown cells 
was assayed by western blot. 

Immuno-electron-microscope and immunostaining 

For EM immunocytochemistry, samples were pre- 
pared according to (37). The cells were fixed in 4% 
formaldehyde/0.2% glutaraldehyde in 0.1 M PHEM 
[60 mM PIPES, 25 mM HEPES, 2 mM MgCl 2 , 10 
mM EGTA (pH 6.9)], cryoprotected in 2.3 M sucrose, 
mounted on aluminum pins and frozen in liquid nitro- 
gen. Thin-frozen sections were then cut on a Leica EM 
UC6/FC6 cryomicrotome (Leica), collected on a drop of 
sucrose/methylcellulose mixture and placed on a formvar- 
carbon grid. The sections were labeled with primary 
antibody and the label was subsequently visualized by 
colloidal gold conjugated to protein A. Sections were 
stained/embedded in 2% methylcellulose/0.2% uranyl 
acetate and observed under a Tecnai 12 transmission 
electron microscope. Immunostaining was performed as 
described (38). 

Expression vector and cell transfection 

Cloning vector pBKSII (BKS) is commercially available 
from Agilent/ Strategene, Inc. Short hairpin RNA (shRNA) 
expression vectors were generated as described previously 
(39). Short hairpins of 27-29 nucleotides for PARP-1 were 
expressed by the U6 promoter in the pENTR/DTopo 
vector (Invitrogen). PGEM1 plasmid containing the U6 
promoter was used as the template for PCR with the 
Sp6 primer (CACCGATTTAGGTGACACTATAG) and 
PARP-1 primer (AAAAAAAG GACAAGAC GTACGC 
TAA GAACAAC CCCCAAG CTTCGGAG TTGTTC 
TTA GCGCACA TCTTGTCC CGGTGT TTCGTCCTT 
TCCACAA) (40). Cells were transfected by Effectene trans- 
fection reagent (Qiagen) according to manufacturer's in- 
structions. 

Chemicals and cell treatment 

To synchronize cells, nocodazole (Sigma) was dissolved 
at 10 mg/ml in dimethyl sulfoxide (DMSO) and stored 
at -20°C. PARP-1 inhibitor, olaparib (AZD-2281, Sell- 
eckchem), was dissolved in DMSO (10 mM stock solution) 
and stored at — 20° C. EGF (Sigma) was dissolved in phos- 
phate buffered saline (PBS) and stored (10 |xg/ml stock) at 
4°C. Small aliquots of these stock solutions were used once. 



Cell synchronization 

The cells were initially plated at a density of 5 x 10 5 cells in 
a 10-cm dish. One day after seeding, the cells were washed 
twice with prewarmed PBS at 37°C and incubated with or 
without 100 ng/ml nocodazole-containing medium. Eigh- 
teen hours later, the media were removed and the cells were 
washed twice with prewarmed PBS. Mitotic arrested cells 
were then harvested and used either for western blotting or 
ChlP. 

Olaparib and EGF treatment 

Nocodazole-treated cells were released by changing fresh 
media after washing twice with PBS in the presence and ab- 
sence of olaparib (10 |xM) for 2 h, unless otherwise indi- 
cated, at 37° C in culture dishes, followed by harvesting the 
cells. For the EGF experiment, cells were treated with EGF 
(50 ng/ml) for 1 h (41). Cells were treated with nocodazole 
or olaparib or EGF in 60 mm dishes. 

Flow cytometric analysis of cell cycle status 

After nocodazole treatment, control or shRNA PARP_1 trans- 
fected cells were released at different time points (0, 1:30, 
2:30, 4 and 5:30 h), then cells were harvested, washed twice 
with ice-cold PBS and fixed in 70% ethanol. Tubes contain- 
ing the cell pellets were stored at — 20° C. Afterward, the 
cells were centrifuged at 1000 g for 10 min and the super- 
natant was discarded. The pellets were resuspended in 30 
(xl of phosphate/citrate buffer (0.2 Na2HPC>4 and 0.1 citric 
acid, pH 7.5) at room temperature for 30 min. Cells were 
then washed with 5 ml of PBS and incubated with 400 |xl 
of propidium iodide (PI) solution (50 |xg/ml PI, 1 mg/ml 
RNase A and 1 g glucose/1 of PBS) for 30 min. The sam- 
ples were analyzed on a Coulter Elite flow cytometer. 

Antibodies 

The following antibodies were used for the ChlP experi- 
ment: rabbit polyclonal anti-PARP (ab6079, Abeam), anti- 
H2AX (#39689, Active motif), anti-H2A.Z (#39943, Active 
motif), mouse anti-NFATC2 (ab2722, Abeam) and anti- 
H2A (Active motif, #39111). For western blot analysis, 
anti-PARP- 1 (C2-10, Trevigen), anti-H3 (Millipore), anti- 
H3S10 (Millipore), anti-H3T3 (Millipore), mouse mon- 
oclonal anti-a tubulin (Sigma), anti-NFATC2 (Thermo, 
#PA5-19164), anti-RNAPolII (Covance, 8WG16) and anti- 
pADPr (Trevigen, 10H) were used. For electron microscopy 
(EM) immunocytochemistry, anti-Hl antibody (sc8030, 
Santa Cruz) and anti-PARP- 1 (ab6079, Abeam) were used. 
Either goat anti-rabbit or anti-mouse secondary antibodies 
conjugated to horseradish peroxidase (Sigma) were used. 

Western blotting 

Western blotting was performed on nocodazole-arrested re- 
leased shRNA contro1 cells or shRNA PARP_1 -transfected cells 
using the detection kit from Amersham/GE Healthcare 
(#RPN2106), according to the manufacturer's instructions. 
For semi-quantitative protein analysis, whole cell extracts 
(0.8 x 10 5 cells) were prepared by boiling cells for 10 min 
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in sodium dodecyl sulfate (SDS) sample buffer [25 mM Tris 
(pH 6.8), 2% p-mercaptoethanol, 3% SDS, 0. 1% bromophe- 
nol blue, 5% glycerol] at 1 x 10 7 cells/ml. Lysed cell extracts 
were resolved by SDS-PAGE (4-10% NuPAGE, Invitrogen) 
and transferred to nitro-cellulose membrane (BioRad) by 
i-Blot dry blotting system (Invitrogen). Detection was per- 
formed with ECL-Plus (Amersham) and autoradiography 
film (HyBlot CL). 

RNA isolation and qPCR 

Total RNA was isolated according to the manufacturer's in- 
structions (Qiagen) from shRNA contro1 and shRNA PARP_1 - 
transfected cells or olaparib-treated cells. Purified total 
RNA was treated by Deoxyribonuclease I (Qiagen). First 
strand of cDNA was synthesized from 5 |xg of purified 
DNAse-treated total RNA, according to manufacturer's in- 
structions (Invitrogen). Quantitative PCR was performed 
on a StepOnePlus q-PCR System (Applied Biosystems) us- 
ing 2x SYBR Green master PCR mix (Applied Biosys- 
tems). All amplifications were performed in triplicate using 
2.0 |xl of complementary DNA (cDNA) per reaction. Trip- 
licate mean values were calculated according to the AACt 
quantification method using GAPDH gene transcription as 
reference for normalization. Changes in expression were 
quantitated by the A A threshold cycle (AACt) method as 
described (42). Control reactions lacking reverse transcrip- 
tase yielded little to no signals. 

RESULTS 

PARP-1 remains bound to chromatin during mitosis 



We first compared the PARP-1 protein distribution in in- 
terphase chromatin and metaphase-arrested chromosomes 
(Figure 1 A), confirming the metaphase status of DNA by 
co-staining with phosphoH3/Serl0, a well-known marker 
of mitotic chromatin (43) (Supplementary Figure SI). 
PARP-1 remained bound to condensed chromatin during 
mitosis (Figure 1A). To assess the distribution PARP-1 in 
metaphase chromatin, we performed ultrastructural anal- 
ysis of metaphase chromosomes using immuno -electron 
microscopy (EM) (37) and immunofluorescence. EM im- 
munocytochemistry deploying anti-PARP-1 and the anti- 
body to linker histone HI (Figures 1B-D) revealed a num- 
ber of well-defined domains in the condensed chromatin 
occupied by PARP-1. Confocal microscopy (Figure IE) 
confirmed that PARP-1 and HI localized to distinct non- 
overlapping blocks in mitotic chromatin (Figure IE). Sim- 
ilar anti-correlation in PARP-1 and HI distributions has 
been previously observed in interphase chromatin (28,35). 

To determine exact genomic distribution of PARP- 1 in 
metaphase chromosomes, we performed ChlP-Seq assays 
with anti-PARP-1 antibody (see 'Materials and Methods' 
section) in both asynchronous and mitosis-arrested cells 
(Supplementary Table SI). We used ChlP-Seq signal with 
Input and IgG as a background control and confirmed the 
specificity of ChlP-Seq protocol by determining PARP-1 
occupancy at previously reported PARP-1 -dependent loci, 
PVALB and GDF15 (35), versus constitutively expressed 
PARP-1 -independent genes, GAPDH and Actin. PARP-1 
binding was enriched at PARP-1 -dependent genes (Figure 
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Figure 1. PARP-1 is associated with chromatin during mitosis. (A) Equal 
amounts of total chromatin proteins extracted from asynchronous cells 
and cells arrested in mitosis were analyzed after PAGE on western blot 
using anti-PARP-1, anti-RNAP (RNA polymerase II), anti-H3 and anti- 
phospho-Serl0-H3 (marker of mitotic chromatin) antibodies (B-D). Panel 
(B) shows typical single mitotic chromosome used for Immuno-EM anal- 
ysis. Immuno-EM analysis of PARP-1 (C) and HI histone (D) localization 
in mitotic chromatin. Scale bars correspond to 1 mm (B) and 100 nm (C 
and D). (E) Confocal microscopy analysis of PARP-1 (red) and HI histone 
(green) localization in mitotic chromatin. Lack of overlap between red and 
green areas indicates the extent to which PARP-1 and HI occupy different 
genomic domains. Scale bars correspond to 5 nm. (F-H) The identification 
of PARP-1 protein binding in mitotic (M) and interphase (asynchronous) 
(A) chromatin using ChlP-Seq assays. UCSC genome browser view is pre- 
sented based on ChlP-Seq data. Red box represents the area of promoter 
sequences. Significant peaks are shown as green. (F) PARP-1 binding in 
the previously reported PARP-1 -dependent gene PVALB. (G) Absence of 
PARP-1 binding along PARP-1 -independent GAPDH locus. (H) PARP-1 
binds along KDM2A gene sequences. 
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Figure 2. Identification of PARP-1 binding sites in mitotic and interphase chromatin in HEK293 cells. (A) PARP-1 binding sites predominantly group 
around transcriptional start sites (TSS) during mitosis (red) genome-wide and spread out more evenly in asynchronous cells (green). (B) PARP-1 binding 
in mitotic cells (red) and asynchronous cells (green). Shared sites in mitotic and asynchronous cells are shown in yellow. (C) Prevalent distance of PARP-1 
binding site overlaps between mitosis and interphase. Red box highlights 200 bp (nucleosomal region) distances. (D) PARP-1 binding is associated with 
non-overlapping sets of genes in mitosis versus interphase: 6048 transcription units have at least one peak in the mitosis stage, 2550 transcription units have 
at least one peak in interphase, and 678 transcription units have peaks in mitosis and interphase (shared). (E) interphase-specific PARP-1 binding correlates 
with gene silencing. Using RNA-seq data, we tested whether the difference between the fraction of PARP-1 -occupied genes among non-expressed genes 
and the fraction of PARP-1 -occupied genes among expressed genes is statistically significant (48). 



1F-H; Supplementary Figure S1E), but completely absent 
at PARP-1 -independent gene sequences in both metaphase 
and interphase chromatin (Figure 1 G; Supplementary Fig- 
ure S1E; Supplementary Figure S2). 



PARP-1 accumulates near transcriptional start sites in mi- 
totic chromatin 

Analysis of cell-phase specific PARP-1 occupancy iden- 
tified three groups of PARP-1 binding sites: interphase- 
specific, mitotic-specific and shared (Figure 2), indicating 
that PARP-1 performs different functions during different 
segments of the cell cycle. We identified 12 823 interphase 
PARP-1 binding sites (Supplementary Table S2) and 25 
388 mitotic sites (Supplementary Table S3) (Supplemen- 
tary Figures S3 and S4). Analysis of PARP-1 occupancy 
relative to annotated genomic features in both mitotic and 
interphase chromatin revealed marked binding enrichment 
near transcriptional start sites (TSSs) (Figure 2A). Based 
on technical limitations of ChlP-seq resolution, any dis- 
tance less than 200 bp was considered below method resolu- 
tion and assumed to represent the same binding site. Only 
135 binding sites that shared between interphase and mi- 



totic chromatin were located closer than 200 bp from one 
another (Supplementary Table S4) (Figure 2B). The over- 
whelming majority of PARP-1 binding sites was different 
between mitotic and interphase chromatin, strongly sug- 
gesting that PARP- 1 becomes redistributed during the tran- 
sition from mitosis to interphase (Figure 2C). 

Next, we examined how PARP-1 occupancy corresponds 
to known genes by comparing its binding sites to annotated 
transcriptional units. While 6048 transcriptional units were 
identified as being occupied by PARP-1 in mitotic chro- 
matin and 2550 were occupied by PARP-1 during inter- 
phase, only 678 transcriptional units were occupied during 
both segments of the cell cycle (Figure 2D; Supplementary 
Table S5). This overlap is significantly less than expected 
by chance (P = 10~ 43 , hypergeometric test), thus indicating 
that PARP-1 is redistributed away from 'mitotic' TSSs to a 
new set of 'interphase' TSSs. These findings strongly sug- 
gest that PARP-1 marks specific TSSs during mitosis and 
that upon completion of mitosis become redistributed to a 
novel non- overlapping set of TSSs. 

The functional significance of PARP-1 redistribution 
from mitosis to interphase is supported by earlier studies 
of the role of PARP-1 in transcription activation. Specif- 
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Figure 3. PARP-1 colocalizes with H2A.Z-bearing nucleosome. (A) 
PARP-1 binding in relation to nucleosome position. 1200 bp DNA se- 
quences centered at PARP-1 peaks were tested to computationally pre- 
dict nucleosome binding affinity (NBA). Average NBA scores for mi- 
totic and interphase cells demonstrate that PARP-1 binding sites coin- 
cide with nucleosome position. (B) PARP-1 binding sites colocalize with 
H2A.Z histone. Three randomly chosen genes from the mitotic-specific 
gene list were subjected to ChIP analysis using anti-H2A.Z, anti-PARP- 
1 and anti-IgG antibody. PARP-1 and H2A.Z are simultaneously present 
at all tested sites. (C) PARP-1 binding site within MED 15 locus colocalizes 
with H2A.Z/H2AX histones, but not with H2A. Arrows indicate positions 
of oligonucleotide primers, which were used in the ChIP assay, along the 
human genome sequence and ChlP-seq PARP-1 peaks (bottom). 



ically, PARP-1 automodifies itself with pADPr during 
PARP-1 -dependent transcriptional activation, leading to 
electrostatic repulsion between the newly assembled 

pADPr and DNA, which results in dissociation of the au- 
tomodified PARP-1 and chromatin (27,30). Consequently, 
continuing PARP-1 presence at TSS regions during inter- 
phase suggests that transcriptional activation has not oc- 
curred and that these loci will remain inactive throughout 
interphase. Indeed, significant number of genes that have 
interphase-specific PARP-1 binding site within ± 1 .5 kb (up- 
stream or downstream) of TSSs were silent (P = 5 x 10 -6 by 
a hypergeometrical test; Figure 2E) in HEK293 transcrip- 
tome (48) demonstrating a remarkable anti-correlation be- 
tween PARP-1 binding in the interphase and gene expres- 
sion (Figure 2E). Therefore, taken together, our results sug- 



gest that the difference in PARP- 1 distribution between mi- 
totic and interphase chromatin is caused by its automodi- 
fication following gene activation after mitosis, while silent 
genes keep PARP-1 bound to promoters. 

PARP-1 colocalizes with H2A.Z-bearing nucleosome in mi- 
totic chromatin 

Computation prediction of nucleosome binding affinities 
in the areas centered at PARP-1 peaks suggests 100% co- 
occurrence between PARP-1 binding and the presence of 
a nucleosome (Figure 3A). Unfortunately, modeling un- 
ambiguous predictions of nucleosome positioning near the 
TSS is problematic due to deviations from regular nucleo- 
somal arrays caused by Pol II and general transcription ma- 
chinery recruitment. Therefore, we determined experimen- 
tally the localizations of PARP-1 and TSS-associated mi- 
crosomes using ChIP with PARP-1 and core histones an- 
tibodies. Earlier research on Drosophila demonstrated that 
mononucleosomes represent binding platforms for PARP- 
1 (49). The nucleosomes containing histone variant H2Av, 
a homologue of mammalian H2A.Z and H2AX, are pre- 
ferred sites for PARP-1 binding in vivo (30,50). Consis- 
tent with these findings, we observed PARP-1/H2A.Z co- 
occupancy at TSSs of all tested mitotic PARP-1 binding 
sites in HEK293 cells (Figure 3B). The profile of PARP-1 
binding corresponds perfectly to the distributions of H2A.Z 
and H2AX histone variants, exclusive of H2A histone lo- 
calization (Figure 3C and Supplementary Figure S5). Taken 
together, these findings suggest that H2A.Z/X nucleosomes 
are colocalized with PARP- 1 in promoter regions in mam- 
malian genome. 

PARP-1 binding sites co-occur with NFATC2 transcriptional 
factor binding sites genome-wide 

To address which specific genes are marked by PARP-1 
during mitosis, we performed consensus sequence compar- 
ison for the mitotic and interphase sites of PARP-1 occu- 
pancy. This analysis revealed GGAAAGG consensus se- 
quences in proximity to approximately 65% of all sites of 
PARP-1 localization in mitosis and interphase (Figure 4A 
and Supplementary Figure S6). This sequence is similar to 
the known binding sites of transcriptional factor NFATC2 
(51). NFATC2 serves as a transcriptional factor control- 
ling inducible expression of interleukin gene family (52,53). 
Recent research has shown that transcription factors from 
NFAT family act as regulators of oncogenic transforma- 
tion (54,55). NFATC2 remains bound to chromatin dur- 
ing interphase, as does PARP-1 (Figure 4B). We deter- 
mined NFATC2 occupancy at a randomly selected group 
of PARP-1 -positive peaks. Indeed, three tested PARP-1- 
positive loci OR51B2, ZSCAN4 and RRAS2 demonstrated 
significant occupancy by NFATC2, while the PARP-1 neg- 
ative locus TMCC2 showed no NFATC2 binding in ei- 
ther mitosis or interphase (Figure 4C-D). Genome-wide 
co-occurrence of PARP-1 and NFATC2 binding sites sug- 
gests that genes marked by PARP-1 may share a com- 
mon mechanism of regulation by NFATC2. This hypoth- 
esis is supported by earlier reports that NFATC2 directly 
interacts with PARP-1 and PARP-1 mediated poly(ADP- 
ribosyl)ation (56,57). 
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Figure 4. PARP-1 binding sites co-occur with NFATC2 binding consensus. (A) Consensus binding motif found in proximity to PARP-1 binding sites. The 
most significant motifs associated with PARP-1 binding sites in mitotic and asynchronized cells were compared to known NFATC2 binding motif. (B) 
NFATC2 remains associated with chromatin during mitosis. Equal amounts of total chromatin proteins extracted from asynchronous and mitosis-arrested 
cells were analyzed by PAGE and western blot using anti-NFATC2, anti-PARP-1, anti-RNAP (RNA polymerase II), anti-H3 and anti-phospho-Serl0-H3 
(marker of mitotic chromatin) antibodies. CF, chromatin fraction; N-CF, non-chromatin fraction. (C) NFATC2 binds to TGGAAA sequence in three 
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Functional classification of genes occupied by PARP-1 

We used Gene Signature Enrichment Analysis (GSEA) to 
identify 33 gene ontology (GO) categories that are enriched 
in at least one PARP-1 occupancy group (see Figure 2; oc- 
cupied in mitosis, but no longer occupied in interphase, 
occupied in both mitosis and interphase, and occupied in 
interphase). We performed GSEA using genes containing 
PARP-1 binding sites and separately with genes that con- 
tained NFATC2 motifs around PARP-1 binding site (Fig- 
ure 5). Clusters of genes that are occupied by PARP-1 
specifically during mitosis (Figure 5, clusters 1, 2 and 3) con- 
tain several enriched GO categories that are related to cell 



adhesion, intercellular interaction, receptor-mediated sig- 
naling and protein secretion. Curiously, very few GO cat- 
egories, including defense response, carbohydrate binding 
and protein kinase A binding, are overrepresented among 
genes that are occupied by PARP- 1 during interphase. Con- 
sistent with functional systems regulated by NFATC2 these 
sites are enriched in GO categories related to the circula- 
tory system (cluster lc), defense response and carbohydrate- 
binding genes (cluster 4). 
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Figure 5. The distribution of PARP- 1 binding sites and NFATC2 consen- 
sus between functional groups of genes in human genome. GO enrichment 
among genes from different categories based on PARP-1 occupancy dur- 
ing mitosis and interphase. Using DAVID software, genes displaying differ- 
ent PARP-1 occupancy patterns during mitosis and interphase were tested 
separately for enrichment with respect to gene ontology terms including 
categories for biological processes (bp), molecular functions (mf) and cel- 
lular components (cc). The occupancy pattern is indicated on the bottom 
as follows: the mitotic group (M) has PARP-1 peak(s) in mitosis within 10 
kb from TSS, mitotic gene in interphase (M gone in I) had PARP-1 peak(s) 
in mitosis within 1.5 kb, but lost them during interphase [this region shown 
to be important for silencing during interphase by PARP-1], mitotic and 
interphase (M and I) has peaks within 10 kb of TSS in both mitosis and 
interphase and interphase (I) has peaks in interphase only. An additional 
GO enrichment test was performed for genes that contained NFATC2 mo- 
tif within 200 bp of the PARP-1 peak center. Intensity of red shading on 
the heatmap corresponds to the degree of significance enrichment. 



Bookmarking promoters with PARP-1 during mitosis is re- 
quired for transcription restart after mitosis completion 

Our data demonstrate that PARP-1 binds promoters of 
a subset of genes during mitosis. Following mitosis, the 
amount of pADPr in chromatin increases dramatically 
(Supplementary Figure S7A). Therefore, it seems likely that 
PARP-1 bookmarks specific genes during mitosis and con- 
tributes their transcriptional reactivation by poly(ADP- 
ribosyl)ation of chromatin post-mitotically. To test whether 
PARP-1 is required for transcriptional restart after mito- 
sis, we compared primary transcript levels of PARP- 1 as- 
sociated genes in PARP-1 knockdown and control cells af- 
ter mitosis. We transfected HEK293 cells with shRNA contro1 
and shRNA PARP_1 using Effectine transfection reagent (Qi- 
agen) for 36 h. To synchronize cells in mitosis we treated 
cells with nocodazole (100 ng/ml) for 18 h. Synchronized 
cells were released by washing with PBS and culturing for 
2 h. The protein levels of PARP-1 were evaluated by west- 



ern blotting. shRNA transfection resulted in 65% PARP- 
1 depletion and reduced total pADPr by 88% (Figure 6A 
and B), but did not affect cell cycle progression (Supple- 
mentary Figure S7B-C). We found that a shorter treatment 
with shRNA failed to sufficiently deplete PARP-1 because 
of high stability of the PARP-1 protein. A prolonged treat- 
ment results in a greater PARP-1 depletion but affected cell 
viability due to secondary effects. 

We compared post-mitotic levels of primary transcripts 
(Figure 6A and B) for two groups of genes selected from 
the pool of PARP-1 mitotic- specific sites: randomly selected 
genes (MED15, LAM CI, ZSCAN4 and KDM2A) (Supple- 
mentary Figure S7D) and EGF hormone-dependent genes 
(FOS, EGR2, LDLR, BRCA1 and EGR3) (41). To distin- 
guish total mRNA from new (nascent) RNA, we designed 
two pairs of primers for each gene: exon to exon (EE) and 
intron to exon (IE) (Figure 6C). Newly transcribed RNA 
(IE) cannot be detected during mitosis. After mitosis, the 
amount of nascent RNA increases rapidly for PARP-1- 
occupied genes, but slowly for a PARP-1 -independent gene 
(GAPDH) (Supplementary Figure S7E). 

Following mitosis, the total amount of mRNA detected 
by qPCR using EE primers remained the same (Figure 6D), 
as was evidenced by the unchanged volume of cytoplasm. 
However, the production of new RNA, as detected by IE 
primers, was reduced by at least two-fold in PARP- 1 knock- 
down cells relative to control cells (Figure 6D,E). A simi- 
lar decline in the level of nascent transcripts was observed 
in cells pretreated with the specific PARP- 1 inhibitor (Ola- 
parib) instead of depleting PARP-1 using shRNA (Supple- 
mentary Figure S8 and 6F-G). Thus, PARP-1 binding to 
TSSs does facilitate transcriptional reactivation of PARP-1 
associated genes after mitosis by a mechanism that relies on 
PARP- 1 -dependent poly(ADP-ribosyl)ation. 

DISCUSSION 

The condensed state of mitotic chromatin poses a chal- 
lenge for transmission of epigenetic information across cell 
division, as most transcriptional factors and RNA poly- 
merases dissociate from mitotic chromatin during this pro- 
cess (2,6). The few factors that remain attached must be suf- 
ficient to restore proper active chromatin patterns upon the 
completion of cell division. Here, we report a new PARP-1- 
dependent bookmarking mechanism of mitotic chromatin. 
Our findings show PARP- 1 acts as an epigenetic marker for 
a large subset of genes. Our data strongly suggest that being 
stably bound to transcription start sites, PARP- 1 can trigger 
quick and robust transcription restart of these genes follow- 
ing the completion of mitosis. We demonstrate that effector 
proteins, such as PARP- 1 can serve as epigenetic marks that 
remain stably associated with certain promoters, along with 
posttranslational modification of histones or DNA. 

Previously, activation of PARP- 1 has been implicated in 
stimulating transcription by loosening chromatin during 
developmental processes and in response to environmental 
stress (26-28,50,58). Upon PARP-1 activation, poly(ADP- 
ribose), a nucleic acid analogue which is twice as nega- 
tively charged as DNA, is synthesized from NAD+ and 
either form long polymers or modifies acceptor proteins 
(59). The negatively charged poly(ADP-ribose) then par- 
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Figure 6. PARP-1 and its enzymatic activity are required to restart post-mitotic transcription. (A and B) Generating synchronous PARP-1 knockdown 
cells (Supplementary Figure S7B). PARP-1 level and enzymatic activity in PARP-1 shRNA and control shRNA-expressing cells. Each experiment was 
repeated three times. (A) Mouse monoclonal anti-PARP-1 and anti-pADPr (10H) antibodies were used to detect the presence of PARP-1 protein and 
poly(ADP) ribose moieties. Anti-tubulin antibody was used as a loading control. Quantification of PARP-1 and pADPr amounts based on three repeated 
experiments was performed using ImageQuant software (B). (C) Two types of primer pairs were employed to detect total mRNA (EE) and nascent RNA 
(IE). Arrows represent primers. (D and E) Comparative analyses of transcriptional re-start after mitosis in control and PARP-1 depleted cells for randomly 
chosen PARP-1 -dependent genes (D) and EGF-PARP-dependent genes (E). Each experiment was repeated three times. The fold difference in q-RT-PCR 
signal between PARP-1 knockdowns and control is marked along the j-axis. Transcription of EGF-dependent genes was stimulated by EGF treatment 
as described elsewhere (41). (F) PARP-1 inhibition in human cells (Supplementary Figure S8). Comparative analysis of PARP-1 activity (pADPr) in cells 
treated with PARP-1 inhibitor, olaparib and in untreated control cells. Mouse monoclonal anti-pADPr (10H) was used to detect pADPr. Anti-tubulin 
antibody was used as a loading control. (G) Comparative analyses of transcriptional re-start after mitosis in cells treated with PARP-1 inhibitor, olaparib 
and in untreated control cells, for randomly chosen PARP-1 -dependent genes. Each experiment was repeated three times. The fold difference in q-RT-PCR 
signal versus input between control and olaparib-treated samples is shown along the 7-axis. 



ticipate in electrorepulsive shuttling of chromatin proteins, 
such as histones and ultimately lead to their dissociation 
from DNA, which results in chromatin loosening and in- 
creased accessibility of DNA for transcription factors and 
general transcription machinery (26,27,30,60,61). 

In agreement with these findings, we observed that mi- 
totic and interphase PARP- 1 cistromes are essentially non- 
overlapping and that PARP-1 becomes redistributed in 
chromatin following completion of mitosis, leaving promot- 
ers of genes that are associated with PARP- 1 driven activa- 
tion. Intriguingly, the retention of PARP-1 near TSSs dur- 
ing interphase strongly correlates with gene silencing. 

The localization and dynamics of PARP- 1 in mitotic and 
interphase chromatin is consistent with PARP-1 acting as 
an epigenetic memory mark. Indeed, (i) PARP-1 remains 
stably bound to TSS of a specific subsets of genes, (ii) 
PARP- 1 correlates with the renewed transcription of these 
genes upon mitosis exit and (iii) PARP-1 dissociates from 
these genes upon restart of transcription. As PARP-1 disso- 



ciates from many of its targets during interphase, the use of 
mitotic chromatin for genome-wide identification of PARP- 
1 -bound promoters produced a more comprehensive list 
of PARP-1 -dependent genes, as compared to screening in- 
terphase chromatin of asynchronous cells. In that respect 
PARP- 1 acts as an active epigenetic mark that can regulate 
its own catalytic activity and the state and accessibility of 
surrounding chromatin (Figure 7), as opposed to passive 
marks such as histone modifications that require recruit- 
ment of additional modification-sensing effectors. 

Remarkably, PARP-1 is required for the expression of 
several genes that are responsible for products that are in- 
volved in interactions between cells and with extracellular 
matrix. As mitotic cells typically lose contact with other 
cells and the extracellular matrix, marking genes respon- 
sible for establishing cellular identity by an epigenetic fac- 
tor, is necessary for rapid restoration of cell niche inter- 
actions, preventing uncontrolled cell migration and poten- 
tially tumorigenic de-differentiation. Upon the completion 
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Figure 7. PARP-1 acts as an active epigenetic memory marker to regulate post-mitotic transcription. Based on our finding we propose the following model 
of PARP-1 action as an epigenetic mark: PARP-1 remains bound to chromatin during mitosis and facilitates transcription after mitosis by poly(ADP- 
ribosyl)ating chromatin, effectively opening it. Upon depletion or inhibition of PARP-1, chromatin remains condensed after mitosis, thereby blocking or 
delaying transcriptional restart. Green triangle indicates NFATC2 transcription factor binding in proximity to PARP-1 binding site. 



of mitosis, PARP-1 opens up chromatin around TSSs of 
these genes, thereby allowing access to transcription ma- 
chinery and reactivating transcription. PARP-1 associa- 
tion with genes involved in metabolic and developmen- 
tal pathways, cell adhesion, communication, proliferation 
and apoptosis suggests that PARP-1 may serve as a central 
developmental/pathogenic hub, controlling such basic or- 
ganismal functions as cellular proliferation and adhesion, 
while deregulation of its production and localization can 
lead to pathogenesis. 
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